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In the present paper we obtain the WIMP velocity distribution in our vicinity starting from 
spherically symmetric WIMP density profiles in a self consistent way by employing the Eddington 
approach. By adding a reasonable angular momentum dependent term in the expression of the 
energy, we obtain axially symmetric WIMP velocity distributions as well. We find that some density 
profiles lead to approximate Maxwell-Boltzmann distributions, which are automatically defined in 
... a finite domain, i.e. the escape velocity need not be put by hand. The role of such distributions 

0^ ' in obtaining the direct WIMP detection rates, including the modulation, is studied in some detail 

, and, in particular, the role of the asymmetry is explored. 
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INTRODUCTION 



The combined MAXIMA-1 [H,[i,[l, BOOMERANG [3,[3, DASI ^ and COBE/DMR Cosmic 
Microwave Background (CMB) observations imply that the Universe is flat fs*], Q = 1.11 ± 0.07 
and that most of the matter in the Universe is Dark Q. i.e. exotic. Combining the recent WMAP 
data llQ] with other experiments one finds: 



nt = 0.0456 ± 0.0015, ilcDM = 0.228 ± 0.013, = 0.726 ± 0.015 



Since the non exotic component cannot exceed 40% of the CDM [ll|, there is room for the exotic 
WIMP's (Weakly Interacting Massive Particles). Supersymmetry naturally provides candidates 
for the dark matter constituents [l^-[l3j|. In the most favored scenario of supersymmetry the 
^vq LSP (Lightest Supersymmetric Particle) can be simply described as a Majorana fermion, a linear 
00 ^ combination of the neutral components of the gauginos and higgsinos [l3| j Hi ' 01 ' [3 > 12 1- most 
' calculations the neutralino is assumed to be primarily a gaugino, usually a bino. Other particle 
', models have also been considered, like Kaluza-Klein WIMPs (see, e.g. the recent work ^§1 and 
references there in ), sterile neutrinos [l3|, technicolor 20] and recently composite WIMP's |2l| (see 
also the recent theory review ^2^). 

Even though there exists firm indirect evidence for a halo of dark matter in galaxies from the 
observed rotational curves, it is essential to directly detect (T^j-ll^l such matter. Until dark matter 
is actually detected, we shall not be able to exclude the possibility that the rotation curves result 
from a modification of the laws of nature as we currently view them. This makes it imperative that 
5J] ' we invest a maximum effort in attempting to detect dark matter whenever it is possible. Furthermore 
5^ ' such a direct detection will also unravel the nature of the constituents of dark matter. 

The possibility of such detection, however, depends on the nature of the dark matter constituents 
(WIMPs). Since the WIMP is expected to be very massive, ruy, > 30GeV, and extremely non 
relativistic with average kinetic energy T « 50KeV{m^/100GeV), it can be directly detected [isj - 
[23] mainly via the recoiling of a nucleus (A,Z) in elastic scattering. The event rate for such a process 
can be computed from the following ingredients: 

1. An effective Lagrangian at the elementary par ticle (qua rk) level obtained in the framework of 



supersymmetry as described , e.g., in Refs |14j-jl7l|.j24|. 
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A well defined procedure for transforming the amplitude obtained using the previous effective 
Lagrangian from the quark to the nucleon level, i.e. a quark model for the nucleon. This step is 
not trivial, since the obtained results depend crucially on the content of the nucleon in quarks 
other than u and d. This is particularly true for the scalar couplings, which are proportional 



to the quark masses |25j — [26| , [27[ as well as the isoscalar axial coupling [27j, |28 



3. Knowledge of the relevant nuclear matrix elements |29| , [30| , [31| , obtained with as reliable as 
possible many body nuclear wave functions. Fortunately in the case of the scalar coupling, 
which is viewed as the most important, the situation is a bit simpler, since then one needs 
only the nuclear form factor. 

4. Knowledge of the WIMP density in our vicinity and its velocity distribution. Since the essential 
input here comes from the rotational curves, dark matter candidates other than the LSP 
(neutralino) are also characterized by similar parameters. 

In the past various velocity distributions have been considered. The one most used is the isother- 
mal Maxwell-Boltzmann velocity distribution with < >= (3/2)ug where vq is the velocity of the 
sun around the galaxy, i.e. 220 km/ s. Extensions of this M-B distribution were also considered, in 
particular those that were axially symmetric with enhanced dispersion in the galactocentric direc- 
tion [32, 3^ l^]. In all such distributions an upper cutoff Vesc = 2.84t'o was introduced by hand. 



in the range obtained by Kochanek[35|. In a different approach Tsallis type functions, derived from 
simulations of dark matter densities were employed, see e.g. recent calculations ,36] and references 
there in . 

Non isothermal models have also been considered. Am ong those one should mention the late infall 
of dark matter into the galaxy, i.e caustic rings [37, H, 33, 3, 41 1, dark matter orbiting the Sun 



[42| . Sagittarius dark matter [43|. 

The correct approach in our view is to consider the Eddington proposal i.e. to obtain both the 
density and the velocity distribution from a mass distribution, which depends both on the velocity 
and the gravitational potential. Our motivation in using Eddington [43| approach to describing the 
density of dark matter is found, of course, in his success in describing the density of stars in globular 
clusters. Since this approach adequately describes the distribution of stars in a globular cluster in 
which the main interaction is gravitational and because of its generality , we see no reason why such 
an approach should not be applicable to dark matter that also interact gravitationally. It seems, 
therefore, not surprising that this approach has been used by Merritt [3| and applied to dark matter 
by UUio and Kamionkowski [iif and more recently by us [i3|) [Sl- 
it is the purpose of the present paper to extend the previous work and obtain a dark matter 
velocity distribution, which need not be spherically symmetric, consistent with assumed halo matter 
distributions with a natural upper velocity cut off. It will then be shown that this distribution can 
be approximated by a Maxwell-Boltzmann (M-B) distribution with a finite domain that depends 
on the asymmetry parameter. The distribution obtained will be used to calculate WIMP direct 
detection rates including the annual modulation, as a function of the asymmetry parameter. 



THE DARK MATTER DISTRIBUTION IN THE CONTEXT OF THE EDDINGTON 

APPROACH 

As we have seen in the introduction the matter distribution can be given [i^ as follows 

dM = 2tt f{^{r),Vr,Vt) dx dy dz vt dvt dvr (1) 

where the function / the distribution function, which depends on r through the potential $(r) and 
the tangential and radial velocities Vt and Vr- We will limit ourselves in spherically symmetric 
systems. Then the density of matter p(|r|) satisfies the equation: 

dp = 2tt /($(|r|),u,.,wt) Vt dvt dvr (2) 

The distribution is a function of the total energy: 



• The energy E is $(r) + ^i-. Then 



p(r) ^ An /($(r) + —)v^dv 



An / f{E)^/2{E ~^)dE 



This is an integral equation ol the Abel type. It can be inverted to yield: 

^(^) d r dp 



An^ dE Je VW^~E d^ 



The above equation can be rewritten as: 
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(3) 



(4) 



(5) 



In order to proceed it is necessary to know the density as a function of the potential. In practice 
only in few cases this can be done analytically. This, however, is not a problem, since this 
function can be given parametrically by the set {p{r), ^{r)) with the position r as a parameter. 
The potential $(r) for a given density p{r) is obtained by solving Poisson's equation. 
Once the function f(E) is known we can obtain the needed velocity distribution fr^ (v) in our 
vicinity (r = r^) by writing: 



/r.(^^) = /(*WU.. + y)Mr = r,) 



(6) 



• We suppose now that there is an additional kinetic term associated with an angular momentum 
[iil, i.e. models of the Opsikov-Merritt type [H^l , [H , [Hlj : 
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where and Vt are the radial, i.e. outwards from the center of the galaxy, and the tangential 
components of the velocity and tq is the "anisotropy radius" to be treated as a phenomeno- 
logical parameter. We now have: 



p^2n J /($(r 
The last integral takes the form: 

p(r) ^An[l 
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1 



Vt dvt dVr 



f{E)^2{E~^)dE 



This equation is formally the same with Eq. ([3]) with the understanding that 



p($)^p(<i>,ro) 1 + 



{r{4>)Y 



(8) 



(9) 



(10) 



The r($) can be obtained by inverting the equation $ = $(r). In practice this is not needed, 
if, as we have mentioned above, we use r as a parameter. We thus find: 



f{E) 



1 



2\/27r2 



d^ (Pp 1 dp I 



^<A> -E 

The velocity distribution in our vicinity becomes 

/.^(z;)=/(ci>(r)U.^+^ + (l + 4)^) 



(11) 



(12) 



and is only axially symmetric. The isotropic case follows as a special case in the limit ro 0. 

The characteristic feature of this approach is that the velocity distribution vanishes outside a given 
region specified by a cut off velocity Vm, by the positive root of the equation fr,iv) = 



SIMPLE REALISTIC DARK MATTER DENSITY PROFILES 



There are many halo density profiles, which have been employed , see e. g. a recent summary 52 1 
and references therein. Among the most commonly used analytic profiles [53| we will consider the 
following: 

• The VO density profile ^ 



X < c 



P{x) = Po 



' 2(c^ + l) _ (f+lf 



X > c 



(13) 



with a the radius of the Galaxy. The distance c is very large, so that the the rotational velocity 
remains essentially constant with the distance from the center of the galaxy even at quite large 
distances. The above form was taken by the requirement that at a; = c the density is continuous 
with a continuous derivative. This will be referred to as VO profile. The resulting potential 
can be obtained by solving Poisson's equation 



Another simple profile is: 



p{x) 



Po 



X(l+X)2' 



(14) 



known as NFW distribution [5j|,[46[, which has been suggested by N-body simulations . This 
profile provides a better description of the expected density near the center of the galaxy. It 
does not, however, predict the constancy of the rotational velocities at large distances. 

As we have seen above the dependence of the density on the potential is much more interesting. 
With thus show these functions in Figs[T] The obtained rotational velocity curve, due to dark mater 





(b) 

FIG. 1: The NFW density ^ as a function of the potential ^ (a). The same quantity in the case 



of the VO profile (b). $0 = 47rGjva^Po 



alone, is given in Fig. [2] in units of with $0 = AttGno^ Po- One can obtain the value of po 
by fitting the rotational velocity in our vicinity to the sun's rotational velocity 220 km s^^. Using 
a = 3.1 X 10^° m and Xs = 0.8, i.e. our location r., — 0.8a, we find 0.3 and 0.5 GeV/c^cm^"^ for the 
VO and NFW profiles respectively compared to the canonical value [H^ of 0.3 GeV/c^cm^''. 



VELOCITY DISTRIBUTIONS 



Before proceeding further we define the asymmetry in the usual way: 




(a) 



(b) 



FIG. 2: The rotational velocity due to dark matter as a function of the distance in units of 
y^AirGNa^Po- Shown on the left is the one obtained with the VO density profile of Eq. (fT3l) , while 
on the right the NFW profile [i^, see Eq. (HU, was employed. 



and it is usually assumed to be positive, < /? < 0.5. The asymmetry in our vicinity r 
related to the parameter rp introduced above via the relation. 



ro=r.^i^ (16) 

Substitute the last expression in /5(<i>, rg) we get p{^,f3). Viewed as a function of the potential 
p($,/3) and its derivatives up to third order vanish at the zero of the potential. 
The velocity distribution can be cast in the form: 

/..(2/,C,/3) = /($(^.) + ^y'^Y^), ^ = (1^) 

The angle 9 here is defined with respect to a polar axis taken in the direction of the radial component 
of the velocity. 

To get an idea of the effect of the asymmetry we integrate the normalized function over the angles. 
We thus get the distributions shown in Fig. [31 We see that the maximum allowed velocity for the 
NFW distribution we find = 5.2 x 10^ km/s, i.e. is v^ax = 1.21Sqrt\^o\ = 1.215.2 x 10^ 

km/s=6.3 X lO'^km/s, which leads to v^sc — 2.86uoj almost the same with the value previously 
used. On the other hand in the case of the VO profile one finds [i^ \/[^oT — 2-67 x 10^ km/s, i.e. 

Vesc = 3.4uo. 

Since, as we have already mentioned, in many dark matter calculations of the event rates the simple 
symmetric or asymmetric M-B velocity distribution has been employed, one would like to know to 
what extend this can be consistent with standard density profiles [56| . In the present work we would 
like, in particular, to know to what extent such an approximate axially symmetric M-B distribution, 
which vanishes outside the range of y given by the exact distribution, can be derived in the context 
of the above density profiles, and can then be used in calculations related to dark matter searches. 
We seek a solution of the form: 

.^ilil^fKiifl!) 

f^^^ = Ce ^^^(^y (18) 

where C is a normalization constant. The above form has been motivated not only by the structure 
of Eq. (|17p . but in addition by our previous work, in which this shape has been obtained by 
an appropriate limit of velocity distributions described by the radial and tangential Tsallis type 
functions [36.]. The width of the M-B distribution depends on the assumed density profile, following 
the trend of the exact distributions. The parameter b and the allowed range of y, which arc functions 
of /3 are shown in table HI 

Integrating the distribution of Eq. (fTS]) over the angles we obtain the results shown in Fig. HI 




FIG. 3: The velocity distribution in our vicinity, iTTy'^ fi,{y, (3) obtained with the NFW profile as a 
function of the velocity in units of v^j^oj = 5.2 x 10'^ km s on the left and the VO profile with 
-\/$o = 2.67 X 10"^ km s on the right. In the plots the thick, dot-dot, long dash, fine continuous, 
short dash and dot-dash lines correspond to /? = 0.50, 0.39, 0.31, 0.22 0.14 and respectively. 



TABLE I: We show the parameters describing the M-B distribution, which is a good fit to the 
velocity distribution derived from the NFW and VO dark matter density profiles via the 
Eddington approach. The range of y given is in units of y^^o]. 
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FIG. 4: The same as in Fig. [3] but using the approximate M-B velocity distribution. 



TRANSFORMATION INTO OUR LOCAL COORDINATES. 

We must now transform the above distributions from the galactic to the local coordinates. 

fr. (y, /?) fl^^^'hy, 0, 0, /?, a, S, 7) (19) 



This is accomplished by the substitutions: 



where 



y^^X^ + Y^ + Z^ ^,^Y,^2, (20) 



1 1 

X = — (y cos sin 9-^-5 sin a), Y = — {y sin 9 sin (f> — S cos a cos 7) , 
sc sc 



Z = — (ycos6i + (5cosasin7+ 1) , y= — , sc ^ V'^M^ f)i ^2.86 (21) 
sc uo "rot(-^s) 

The angles 6* and (j) are the spherical coordinates defined in the usual way. Here the polar axis has 
been chosen along the sun's direction of motion, the x-axes radially out of the galaxy and the y-axis 
perpendicular in the galactic plane {y — z x x). 6 — 0.135 is the earth's rotational velocity in units 
of the sun's velocity, 7 « 7r/6 is the angle between the axis of the galaxy and the axis of the ecliptic 
and a is the phase of the earth (a = around June the 3nd). 



A BRIEF DISCUSSION OF DIRECT WIMP EVENT RATES (DWER) 



Even though the expressions for the event rates for WIMP detection are well known for the reader's 
convenience we will include the basic formulas here in our own notation [3^, [H^, [H^ : 



(f-) - ^^yP) / J^/i«-l(^,^,0,/^,a,^,7)^c^^- 



The differential cross section is given by 27| : 

du 



da{u, v) = 



[(EsF(w)' + S,p»„Fii(u)] 



where u the energy transfer Q in dimensionless units given by 

M = , Qo = [rupAb]-^ = 40 A--*/^ MeV 
Qo 



(22) 



(23) 



(24) 



with b is the nuclear (harmonic oscillator) size parameter. F(u) is the nuclear form factor and Fii{u) 
is the spin response function associated with the isovector channel. 
The scalar and spin cross sections are given by: 



spin > 



(25) 



where A is the nuclear mass number and /ir(A) (Mr(p) ) is the WIMP-Nucleus (WIMP-nucleon) 
reduced mass. (T ^ ^0 and ct^^q are respectively the WlJVIP-nucleon scalar and spin cross sections 
and C,spin the nuclear spin ME. 

Integrating over the energy transfer u we obtain the event rate for WIMP-nucleus elastic scattering, 
which is given by ^,^,^7, 58]: 



R 



p{0) m 



fcohiA, Hr{A))ap ,^o + fspiniA, firiA))ap'^^o Cspi 
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(26) 



(27) 



IJ-r{p) 



^ *spin (l + ^spincosa) ^^^^ 



The nucleon coherent and spin cross sections are the most important particle physics parameters 
and Cspin is the most important nuclear physics parameter. In this work, however, we are mainly 
interested in the effects of the velocity distribution on the event rate. So to spare the reader the 
inconvenience of detailed discussions of particle and nuclear physics wc will not be concerned with 
the spin cross section. Thus we essentially only need the information of the velocity distribution 
and the effect of the nuclear form factor. Thus the relevant parameters and tcoh, which deal with 
the time averaged rate, and hcoh, which deals with the modulation duo to the annual motion of the 
Earth. They result after the folding of the nuclear form factor with the WIMP velocity distribution. 
More specifically: 



, / dtcoh 11 ^ f dhcoh 1 

tcoh = I — 1 — du , hcoh = - — / —5 — du. 

Jumin du tcoh Jumin dU 



(29) 



The energy transfer u is limited from below by the detector energy threshold and from above by the 
maximum WIMP velocity, i.e. : 

■Wmm = <u< Umax = f , o = v^Hrbvo , Vq = the suu's vclocity (30) 

C^o V a / L J 



One can show that 



dt _ 2 
~ V 3 



^a^F^{u)'^oiaV^), ^ = ^a'^F^{u)H {aV^) cos a. (31) 



where the two functions ^'o(a-\/u) H(ay/u) are respectively the n = (a independent) and the 
cos a Fourier coefficients of the integral: 

J= ydy fr%, 0, cP, f3, a, 6, j)dy (32) 

The integral J contains all the relevant information on the velocity distribution. We see from 
this that essentially the first moment of the velocity distribution appears. Since the high energy 
transfers are suppressed by the nuclear form factor F{u), the above integral is not much affected by 
the behavior of the distribution near the escape velocity. So we expect the M-B approximation to 
be an adequate description of the exact distribution. 

The above expressions manifestly show the essential parts: 

• The elementary cross section. 

This is the most important part. It depends on two ingredients i) The particle model which 
provides the amplitude at the quark level and ii) The procedure for going from the quark to 
the nucleon level. We will not concern ourselves with such issues here. 

• The WIAfP density in our vicinity po- 

This has been obtained in two phenomenological density profile as described above. 

• The dependence of the rate on the cross properties of the target and the WIMP mass. 

• The quantity t. 

This is independent of the parameters of the particle model, except for the WIMP mass. It 
takes into account: 

1. The nuclear structure effects (for the coherent process the nuclear form factor). 

2. The WIMP velocity distribution. 

This may be obtained from a given spherical density profile via the Eddington approach 
as discussed above. 



3. The energy threshold imposed by the detector. In the case of a non zero threshold the 
obtained rates depend on quenching. Such an effect will not be discussed here. The 
interested reader is referred to the literature (see e.g. [59|-[60l|) 

• The quantity h. 

This describes the modulation of the amplitude due to the Earth's motion around the sun. It 
depends on the same parameters as t. 

The evaluation of the quantities t and h with the obtained asymmetric velocity distribution will be 
discussed below. 

With the above ingredients the number of events in time t due to the scalar interaction, which 
leads to coherence 611, can be cast in the form 



R 



-3 



t p(0) m 
ly 0.3GeVcm"3 1kg 280km s"^ 10^6 



s 



-Icoh{A,^lM)) (33) 



where the elementary cross section (7^ ^ocan be treated as a phenomenological parameter. 



RESULTS 



Since the NFW profile preceded the VO profile and is more widely known, we are going to employ 
the velocity distribution obtained with this profile to compute the direct WIMP detection rates. 
We begin with the quantities t and /i, which are pretty independent of the particle model, but are 
sensitive to the WIMP velocity distribution and mass. We first consider differential quantities dt/du, 
which is proportional to the time average differential rate dR/du, and dh/du, which is the relative 
differential modulated rate, as functions of the energy transfer Q. For a typical light and a typical 
heavy target are shown in Figs [5][8l Two values of the WIMP mass, namely 30 and 100 GeV/c^ 
were considered. Note that in some instances the differential modulated rate changes sign. This 
may result in cancellations in the integrated rate yielding smaller values than expected due to the 
smallness of d alone. It may also cause a shift the maximum of the rate from June to December. 




Q — > keV 

(a) (b) 

FIG. 5: On the left we show the differential rate ^^g^ for — 30 GeV in the case of the axially 
symmetric distribution obtained in this work for a heavy nucleus. Here we have taken A=127. On 
the right wc show the same quantity for — 100 GcV. Otherwise the notation is the same as in 
Fig. O The effect of the asymmetry is not visible. For the spin contribution the behavior is 

analogous. 



By integrating the differential rate we obtain the total rates t and /i as a function of the WIMP mass 
for a given energy threshold. Typical examples for the coherent process are shown in Figs [MT21 for 



Q — > keV 

(a) (b) 
FIG. 6: The same as in Fig. [5] for the quantity 




FIG. 7: On the left we show the quantity 



dtr 



for 



30 GeV in the case of the axially 



symmetric distribution obtained in this work for a hght nucleus like A=19. On the right we show 
the same quantity for = 100 GeV. Otherwise the notation is the same as in Fig. [31 For the 

spin contribution the behavior is analogous. 




(a) 



g — > keV 



FIG. 8: The same as in Fig. [7] for the quantity 



du 



(b) 



a typical light target (^^F) and an intermediate-heavy target (^^^I). Two typical threshold values of 
zero and 10 keV were considered. 




— > GeV 

(a) (b) 

FIG. 9: On the left we show the rate tcoh as a function of the WIMP mass in the case of the 
axially symmetric distribution obtained in this work for a heavy nucleus. Here we have taken 
A=127 and assumed zero threshold. On the right we show the same quantity for a threshold 
energy of 5 keV, Otherwise the notation is the same as in Fig. O 




(a) (b) 
FIG. 10: The same as in Fig. [5] for the quantity hcoh- 
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nix — > GeV 

(a) (b) 

FIG. 11: On the left we show the quantity tcoh as a function of the WIMP mass in GeV in the case 
of the axially symmetric distribution obtained in this work for a light nucleus like A=19. On the 
right we show the same quantity for a threshold energy of 5 keV. Otherwise the notation is the 

same as in Fig. |3l 
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FIG. 13: The (time averaged) event rate per kg target per year as a function of the WIMP mass in 
GeV assuming a nucleon cross section = 10~^ pb. The results shown are due to the coherent 
mode for a heavy target (A=127) in (a) and a Ught target in (b) assuming zero threshold . The 
effect of the asymmetry is not clearly visible. At low WIMP masses the rate is suppressed due to 

small reduced mass. 



DISCUSSION 

With the above ingredients we can now compute the total event rate. We do not know what 
the elementary nucleon cross section is, but following the more or less standard practice, we will 
assume it to be independent of the WIMP mass and equal to 10~^pb. We will also take the WIMP 
density in our vicinity to have the canonical value [55|] of 0.3 GeV/c^ cm""^. Employing Eq. ([33| and 
using the values of t discussed above we find the results shown in Fig. [13] for zero energy threshold. 
For an energy threshold value of 10 keV we obtain the results shown in Figs [HI and [T51 for A=127 
and A=19 respectively. On these figures we also show the effect of quenching ^5^]-^] assuming 
an energy threshold of 10 keV. There is no need to show again the effects of modulation on the 
rate, since the above shown results are adequate (h refers to the ratio of the modulated to the time 
averaged rate). 




rrix — > GeV 

(a) (b) 

FIG. 14: The total rates in the case of a heavy target (A=127) computed assuming an energy 
threshold of 10 keV without quenching (a) and with quenching (b) and a nucleon cross section 
ctat — 10^^ pb. Note that the quenched rate is about a factor of 8 smaller compared to that at zero 

threshold (see Fig. [T5)) . 




CONCLUSIONS 



In the present paper we first derived the WIMP velocity distribution for a spherically symmetric 
WIMP density profile. This was done in a self consistent way by applying Eddington's approach. 
Asymmetry of the velocity distribution was also taken into account by incorporating angular mo- 
mentum into the expression of the energy entering the phase space distribution function. Using this 
distribution function we obtained both the differential and total rates entering direct WIMP detec- 
tion. We find that the time averaged differential rates do not sensitively depend on the asymmetry 
parameter (3. The obtained shape of this signal unfortunately cannot really be differentiated from 
that expected from most backgrounds. The differential modulation rate H shows some dependence 
on the asymmetry parameter, especially for light targets. We found that this differential rate changes 
sign at some energy transfer, which depends on the WIMP mass. Perhaps this signature may aid in 
discarding some season dependent backgrounds. 

After that we computed the corresponding total rates. We find that the time averaged total rates 
are not very sensitive to the asymmetry parameter, in agreement with earlier results obtained with 
distributions given in terms of Tsallis type functions [3q |. The maximum rate is expected at a 
WIMP mass of about 75 GeV for a heavy target and 30 GeV for a light target. The modulation 



for a light system is always positive (maximum in June) and tends to increase with asymmetry. For 
zero threshold is rising from 2h = A % to about 2/i = 6 % as the asymmetry increases from /3 = to 
P — 0.5. Assuming a threshold of 5 keV we get 2h — 5 — 8 %. Higher modulations up to 16% can be 
expected at smaller WIMP masses. For a heavy target the modulation increases with asymmetry, 
but it goes through zero at a WMIP mass of about 50 GeV. We expect a positive modulation at 
lower WIMP masses, with a maximum of about 2h = 6% and negative modulation in the case of 
heavier WIMPS ranging to 2|/i| = 4 to 6% depending on the asymmetry parameter. 

Finally we have seen that both the spherically symmetric and the axially symmetric velocity 
distributions obtained from the realistic NFW density profile can be approximated by a Maxwell- 
Boltzmann velocity distribution in a finite domain, i.e. one in which the upper bound of the velocity 
(escape velocity) is not put in by hand but it comes naturally from the Eddington method. Inves- 
tigations of the goodness of such a fit still continue. Anyway, as we have discussed in the section 
marked by DWER on its title, we expect the approximate solution to be appropriate in calculations 
relevant to dark matter searches. So one may use this distribution in the future to simplify the 
calculations. 

We did not consider in this work directional experiments, i.e. experiments in which not only the 



energy but the direction of the recoiling nucleus is also observed |62l | , 63[ . In such experiments both 
the observed time averaged rates as well as the modulation are expected to be direction dependent 
[b^ . We expect that such experiments are going to be much more sensitive to the form of the 
velocity distribution and, in particular, the asymmetry parameter. 
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